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ABSTRACT 

We present results from Suzaku and Swift observations of the nearby radio galaxy 3C 33, and investigate 
the nature of absorption, reflection, and jet production in this source. We model the 0.5-100 keV nuclear 
continuum with a power law that is transmitted either through one or more layers of pc-scale neutral material, 
or through a modestly ionized pc-scale obscurer. The standard signatures of reflection from a neutral accretion 
disk are absent in 3C 33: there is no evidence of a relativistically blurred Fe Ka emission line, and no Compton 
reflection hump above 10 keV. We find the upper limit to the neutral reflection fraction is R < 0.41 for an 
e-folding energy of 1 GeV. We observe a narrow, neutral Fe Ka line, which is likely to originate at least 
2,000 R s from the black hole. We show that the weakness of reflection features in 3C 33 is consistent with 
two interpretations: either the inner accretion flow is highly ionized, or the black-hole spin configuration is 
retrograde with respect to the accreting material. 

Subject headings: galaxies: active - galaxies: jets - galaxies: individual (3C 33) - X-rays: galaxies 



1. INTRODUCTION 

The origin of jets in active galactic nuclei (AGN) is one 
of the most important unsolved problems in extragalactic as- 
trophysics. While 90% of all AGN (Seyfert galaxies and 
radio-quiet quasars) show little or no jet emission, the re- 
maining 10% (the radio-loud AGN and radio-loud quasars) 
launch powerful, relativistic twin jets of particles from their 
cores. Since jets transport a significant fraction of the mass- 
energy liberated during the accretion process, sometimes out 
to ~Mpc distances, understanding how they are produced is 
key to a complete picture of accretion and feedback in AGN. 

X-ray observations of the nuclei of radio-loud and radio- 
quiet AGN are essential for establishing the connection be- 
tween the accretion flow, black hole, and jet. In standard mod- 
els of X-ray emission in AGN, quasi-isotropic power-law X- 
ray emission is generated in a corona that lies above a high 
M ac cretion disk dKrolik & Kallmanl 119871 : iGhisellini etal] 
1994). A fraction of this hard X -ray emission is reflected by 
cold gas in the accretion disk, producing an observed spec- 
trum imprinted with the features of photoelectric absorption, 
fluorescent 6.4 keV Fe Ka line emission, and a Compton 
backscattered reflection conti nuum that produces a bu mp in 
the spectrum above 10 keV dGeor ge & Fab ianl 119911) . The 
Fe Ka line may be relativistically broadened, which is the 
signature of disk reflection near the innermost stable c ircular 
orbit (ISCO) of the black hole (e.g jFabian eTail[l989l) . In ra- 
dio galaxies, there exists an additional power-law component 
of X- ray emission that is associated with the parsec-scale jet 
(e.g. JEvans et~al]|2006l: iBalmaverde et al.ll2006h 

X-ray observations reveal systematic differences between 
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the spectra of radio-loud and radio-quiet AGN. Disk 
reflection is commonly detected in rad io-quiet galaxies 
(Reeves & Turner 2000; R eeves et al.ll2006l) : Compton reflec- 
tion bumps are prominent, and over 50% of local Seyfert 
galaxies show broadened Fe Ka emission from within 25 
Schwarzschild radii. This result may hold even when the 
effects of complex X-ray absorption are taken into ac- 
count dNandra et al.l 120071) . On the other hand, radio- 
loud AGN tend to show weak, or absent neutral Comp- 
ton reflection humps dGrandi et alj|2006t iLarsson et aLll2008t 
S amb runa et al.l2009l : lKataoka et al.l2007l). and in most c ases, 
unresolved neutral Fe Ka lines ((Balmaverde et al. 200^). 

Continuum X-ray observations of radio-loud AGN have 
mostly been restricted to bright broad-line radio galaxies 
(BLRGs) and quasars, which are oriented relatively close 
to the line of sight with respect to the observer. In these 
sources, unabsorbed non-thermal emission from the jet could 
potentially contaminate the unabsorbed accretion-related X- 
ray spectrum and thus dilute the apparent strength of the 
Compton reflection component. Narrow-line radio galaxies 
(NLRGs), on the other hand, which are oriented at low to 
intermediate angles, have the distinct advantage that (unab- 
sorbed) jet-related X-ray emission can be readily spectrally 
separated from (heavily absorbed) accretion-related emission, 
allowing a direct measurement of the strength of Compton re- 
flection. However, narrow-line radio galaxies tend to be rel- 
atively faint X-ray sources compared to broad-line objects, 
meaning that Suzaku observations of these sources are par- 
ticularly useful, owing to the high effective area of the X-ray 
Imaging Spectrometer (XIS) and Hard X-ray Detector (HXD) 
instruments. 

Here, we present the results from Suzaku and Swift obser- 
vations of one of the brightest NLRGs, 3C 33,. It is the only 
NLRG to show potential evid ence for Compton reflection in 
its 2-10 keV X-ray spectrum (Kraft et al. 2007). The source 
may also show evi dence of photoion ized emission lines in 
its X-ray spectrum dTorresi et a l. 2009), which would make it 
one of the few radio-loud AGN to do so. In this paper, we an- 
alyze the Suzaku and Swift spectra (Sj3), and show that reflec- 
tion is not required in the source once the >10 keV spectrum 
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is taken into account (©. We perform detailed diagnostics 
of the Fe Ka bandpass ([Q, and demonstrate that there is no 
evidence for an appreciably velocity broadened component to 
the line. We combine data from non-contemporaneous obser- 
vations with Suzaku, Swift, Chandra, and XMM-Newton, in 
order to search for variability in the source (Sj6]l. We also use 
these data to search for emission lines in the soft X-ray spec- 
trum (§7). We end with a detailed interpretation of the prop- 
erties of absorption, reflection, and the disk-jet connection in 
3C33(®. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Suzaku 

We observed 3C33 with Suzaku on 2007 December 26 
(OBSID 702059010) for a nominal exposure of 100 ks. Both 
the X-ray Imaging Spectrometer (XIS) and Hard X-ray Detec- 
tor (HXD) were operated in their normal modes. The source 
was positioned at the nominal aimpoint of the HXD instru- 
ment. The data were processed using v. 2.1.16 of the Suzaku 
processing pipeline, which includes the latest Charge Trans- 
fer Inefficiency (CTI) correction applied for the XIS. We used 
the standard cleaned events files, which are screened to re- 
move periods during which the satellite passed through the 
South Atlantic Anomaly (SAA), had a pointing direction <5° 
above the Earth, or had Earth day-time elevation angles <20°. 
We describe in detail our analysis of the XIS and HXD data 
below. 

2.1.1. XIS 

We used data from the two operational front-illuminated 
(FT) CCDs (XIS0 and XIS3), together with the back- 
illuminated XIS1 detector. The three XIS CCDs were op- 
erated with a frame time of 8 s. For our analysis, we used 
data taken in the 3x3 and 5x5 edit modes. We selected only 
events corresponding to grades 0, 2, 3, 4, and 6, and removed 
hot and flickering pixels with the CLEANSIS tool. 

We extracted the spectrum of 3C 33 from the XIS CCDs us- 
ing a source-centered circle of radius 2.5', with background 
sampled from an adjacent region free from any unrelated 
sources, as well as the 55 Fe calibration sources at the corners 
of each detector. We generated response matrix files (RMFs) 
for each detector using v. 2007-05-14 of the XISRMFGEN soft- 
ware, and ancillary response files (ARFs) using v. 2008-04-05 
of the XISSIMARFGEN software. 

The net exposure times and count rates for the three CCDs 
are shown in Table Q] We co-added the two FI spectra us- 
ing the ADDASCASPEC program, and grouped the resulting 
spectrum and that of the XIS 1 detector to a minimum of 200 
counts per bin in order to use x 2 statistics. We restricted the 
energy range for our spectral fitting to 0.5-10 keV. 

2.1.2. HXD 

We extracted the source spectrum from the HXD/PIN de- 
tector, using the cleaned PIN events files described above. 
The source was not detected with the GSO instrument. The 
PIN non X-ray background (NXB) spectrum was generated 
from the latest 'tuned' time-dependent instrumental back- 
ground event file provided by the Suzaku Guest Observer Fa- 
cility. We extracted the source and background spectra using 
a common Good Time Interval (GTI) criterion. The source 
spectrum was also dead-time corrected using the HXDDTCOR 
tool. 
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Fig . 1 .— Suzaku HXD/PIN spectrum of 3C 33 (black) and total CXB+NXB 
background spectrum (red). The source spectrum has been binned to 3<r 
above the background level. 




Energy (keV) 

FIG. 2.— Ratio of Suzaku HXD/PIN spectrum of 3C33 to the total 
CXB+NXB background spectrum. Both spectra have been binned using the 
3cr criterion described in the text. 



We estimated the contribution from the Cosmic X-ray 
Background (CXB) using the epoch-dependent PIN re- 
sponse for a flat emission distribution. We used a typi- 
cal C XB spectrum found using HEAO-1 observations (Boldt 
[l987h of the form 9.412 x 10~ 3 x (£7/1 keV)" 129 x 
cxp (—£7/40 keV) photons cm~ 2 s _1 FOV -1 keV -1 , where 
£7 is the photon energy. We added the NXB and CXB 
spectra together using the MATHPHA tool. Using the ap- 
propriate response file for epoch 4 of Suzaku observations, 
ae_hxd_pinhxnome4_20080129.rsp, we binned the source 
spectrum to 3a above the total (NXB + CXB) background 
level. The binned source spectrum and total background spec- 
trum are shown in Figure Q] In Figure [2] we plot the ratio of 
the source spectrum to the NXB + CXB spectrum, using our 
3cr binning criterion. It can be seen that the last bin is only 
marginally significant. We therefore restrict our subsequent 
spectral analysis to 15-30 keV. 

2.2. Swift BAT 

3C 33 is detected in the Swift BAT 22-month all-sky survey 
dTueller et a l. 2009). We used the preprocessed, background- 
subtracted BAT spectrum and standard diagonal response ma- 
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trix made available to the communitjQ The uncertainties are 
calculated from the RMS noise, and include an additional 
multiplicat ive factor to accoun t for systematic noise, as de- 
scribed by lTueller etail (l2009h . Only the first 6 of the 8 total 
channels are above the background, so we discarded the final 
two bins, leaving an energy range for our spectral fitting of 
14-100 keV. 

2.3. Chandra 

3C33 was observed for 20 ks (PI S. S. Murray) with the 
Chandra ACIS CCD camera on 2005 November 08 (OBSID 
6910), and for another 20 ks on 2005 November 12 (OBSID 
7200). Both observations were performed in FAINT mode, 
with the source placed near the standard aimpoint of the S3 
chip. We reprocessed the data using CI AO v4. 1.2 with the 
CALDB v4. 1 .4 to create a new level-2 events file filtered for 
the grades 0, 2, 3, 4, and 6 and with the 0.5" pixel randomiza- 
tion removed. To check for periods of high background, we 
extracted lightcurves for the entire S3 chip, excluding point 
sources. No periods of high background were found in either 
data set. The observation log is shown in Table Q] 

We extracted nuclear spectra from both Chandra data sets 
using the CIAO routine PSEXTRACT. We chose a source- 
centered region of radius 5", and sampled background from a 
surrounding annulus of inner radius 5" and outer radius 10". 
Finally, we grouped the spectra to 25 counts per bin. 

2.4. XMM-Newton 

XMM-Newton observed 3C 33 on 2004 January 21 with the 
EPIC CCDs for 9 ks (PI M. J. Hardcastle). We reprocessed the 
XMM-Newton data using SAS version 5.4.1, and generated 
calibrated event files using the EMCHAIN and EPCHAIN 
scripts. We imposed the additional filtering criteria of se- 
lecting events with only the PATTERN <4 and FLAG=0 at- 
tributes. We searched for periods of high particle background 
by extracting light curves from the whole field of view, ex- 
cluding a circle centered on the source, and selected only 
events with PATTERN=0 and FLAG=0 attributes and over an 
energy range of 10-12 keV (MOS) cameras and 12-14 keV 
(pn). No intervals of high background were found. 

For our spectral analysis we used data from the EPIC pn 
camera only, owing to its greater effective area with respect 
to the MOS cameras. We extracted the source spectrum from 
circle of radius 35" and extracted the background spectrum 
from a large off-source region on the same chip. We grouped 
the resulting spectrum to 30 counts per bin. Finally, we gen- 
erated corresponding RMF and ARF files using the RMFGEN 
and ARFGEN tools. 

3. SPECTRAL ANALYSIS 
3.1. Overview and Motivation of Models 

In this section, we describe our spectral fitting to the broad- 
band X-ray spectrum of 3C33. Since the main goal of this 
paper is to constrain the properties of Compton reflection in 
3C 33, we only consider the Suzaku and Swift data sets here. 
We introduce the Chandra and XMM-Newton spectra in ij6l 
when we perform time-variability analysis of the continuum 
and Fe Ka line. 

We b riefly highlight our spectral fitting methods here. In 
Section l3~2l we show that models with single zones of neu- 
tral absorption (Models I and II) fail to provide a good fit 

6 http://swift.gsfc.nasa.gov/docs/swift/results/bs22mon/ 
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FIG. 3.— Suzaku XIS FI (black), XIS BI (red), HXD PIN (green), and 
Swift BAT (blue) spectra and residuals of 3C33. The model fit is the sum 
of a heavily absorbed power law, a soft, unabsorbed power law (to take into 
account the soft excess), and a Gaussian Fe Ka emission line (Model I). The 
noticeable residuals at ~3 keV motivates us to consider more complicated 
spectral models. 

to the spectrum of 3C33. This motivates us to introduce 
more complex spectral models, such as multiple zo nes o f neu- 
tral absorption (Models III and IV, described in Sj33J and a 
partially io nize d absorbing medium (Models V and VI, de- 
scribed in £|3.4l >. Once a baseline model has been established, 
we attempt to constrain the properties of Compton reflection 
(Models VII and VIII, £j4j>. The casual reader may simply 
wish to consult Table |3]for a summary of the best-fitting mod- 
els. 

In all spectral fits, we linked the parameters of the model 
components across the Suzaku XIS FI, XIS BI, HXD PIN, 
and Swift BAT spectra. We tied the normalizations of the 
FI and BI XIS spectra together, but applied a constant cross- 
normalization factor of 1 .09 for the PIN spectrum relative to 
the XIS as described in the Suzaku Data Reduction Guid^TJ. 
We also included a Gaussian of initial centroid 6.4 keV to 
represent the Fe Ka emission, and allowed its energy, line 
width and normalization to vary freely. We report in detail the 
properties of the Fe Ka line in |5] 

We performed our spectral fitting using v. 1.5. 0-1 3 of the 
ISIS spectral fitting packagfl which has been modified to 
run on a highly parallelized compute platform. All results 
presented here use a cosmology in which tt m ,a = 0.3, Qa.o 
= 0.7, and H = 70 km s" 1 Mpc" 1 . At the redshift of 
3C33, z=0.0597, the luminosity distance is 267.3 Mpc. Er- 
rors quoted are 90 per cent confidence for one parameter of 
interest (i.e., Xmin + 2-7), unless otherwise stated. All spec- 
tral fits include the Galactic absorption to 3C 33 of -/Vh Gal — 
3.03 x 10 20 cm" 2 dKalberla et alJl2005h . 

3.2. Single Zones of Neutral Absorption: Models I and II 

We first fitted to 3C 33 the canonical model for the X-ray 
spectrum of NLRGs: the combination of a heavil y absorbed 
power law and a soft, u nabsorbed power law dEvans et~aT] 
(2006): lBalmaverde et al](l2006h ) (Model I). We tied together 
the photon indices of both power laws. The model resulted 
in a poor fit to the spectrum (x 2 = 242 for 152 dof), with 
noticeable residuals at energies ~3 keV (see Fig. |3J, and an 
atypically low photon index (r=1.39 ± 0.05). An identical fit 

7 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

8 http://space.mit.edu/cxc/isis/ 
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FIG. 4. — As Fig. [5] but fitted to a double partially covered primary law 
(Model III). This model provides a good fit to the spectrum (\ 2 = 176 for 
150 dof). 

is, of course, achieved if we allow the power law to be par- 
tially covered by neutral material (Model II); in this case, the 
covering fraction is 96 ± 1%. 

3.3. Multiple Zones of Neutral Absorption: Models III and IV 

We next investigated the possibility that an additional 
layer of cold absorption is required to fit the broadband 
spectrum of 3C33 adequately. We adopted a double 
partial-covering model, implemented in XSPEC/ISIS as 
phabsxzpcfabs(l)xzpcfabs(2) x (powerlaw+zgauss), where 
ZPCFABS represents the partial coverer (Model III). The in- 
dividual column densities A^h,i and A^h,2, respectively have 
covering fractions fx, f%. 

Table [3] shows the best-fitting parameters for this model, 
and Figure H] shows the counts spectrum, model, and residu- 
als. In short, the model provides a good fit to the spectrum 
(% 2 =176 for 150 dof), and represents a substantial improve- 
ment over Models I and II. Further, the photon index of the 
power law (Ti = 1.71 ± 0.09) is significantly steeper than in 
the previous model, and is now consistent with typical values 
found in narrow-line radio galaxies (e.g., iBalmaver de et alJ 
2006). A gain, an identical fit is found if we choose not to ex- 
plicitly associate the heavily absorbed and unabsorbed spec- 
tral components (Model IV). 

3.4. Warm Absorption: Models V and VI 

An alternative method of modeling the observed spectral 
complexity in 3C 33 is to allow the absorber to be partially 
ionized: such an absorber transmits more of the low-energy 
continuum emission than does a neutral obscurer of identical 
column density. We modeled the Suzaku spectrum as a power 
law, partially covered by a single zone of partially photoion- 
ized absorption (Model V). The ionized absorbed is modeled 
by the ZxiPCF model. 

This model provided a good fit to the spectrum (x 2 = 184 
for 151 dof), with a column density TVh = (5.2t°'^) x 
10 23 cm~ 2 of modestly ionized gas [log(£)=1.5lto2o 
ergs cm s _1 ] as its best-fitting absorption parameters. The 
spectrum and model are shown in Figure Again, an iden- 
tical fit is found if we choose not to explicitly associate 
the heavily absorbed and unabsorbed spectral components 
(Model VI). 

3.5. Definition of baseline model 
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FIG. 5. — As Fig. [5] but fitted to a power law partially covered by a single 
zone of warm absorption (Model V). The model provides a good fit to the 
spectrum (\ 2 = 184 for 151 dof). 

We have demonstrated that the broadband X-ray spectrum 
of 3C 33 can be described with a power law that is partially 
covered by multiple zones of neutral absorption (Model III), 
or by a power law that is modified by a partially-covering 
warm absorber (Model V). Both fits include neutral Fe Ka 
emission. We have also shown that mathematically identical 
fits are achieved if we do not associate the hard and soft X- 
ray components with the same physical process (Model IV 
and Model VI). At this stage, we define as our baseline best- 
fitting spectral models a Model IV, since it is the most statis- 
tically acceptable. We emphasize that this does not change 
the results of our spectral fitting in subsequent sections. We 
provide a detailed discussion of the merits of both the neutral 
and ionized absorption models in [JS] 

4. CONSTRAINTS ON REFLECTION 

We determined the strength of any Compton reflection in 
the spectrum of 3C 33 by replacing the primary power law in 
our baseline model with a PEXRAV component (Model VII). 
This model calculates the sum of the direct power-law emis- 
sion plus the fraction, R, reflected from an infinite slab of 
neutral material. We initially froze the e-folding energy of 
the PEXRAV component at 1000 keV (i.e., we chose there to 
be no spectral cutoff in our Suzaku spectrum), adopted an in- 
clination angle of 60°, and fixed the elemental abundances at 
their solar values. 

This neutral reflection model provided no significant im- 
provement in the fit (x 2 =175 for 149 dof). The best-fitting 
value of the reflection fraction is R = 0, with a 90% confi- 
dence upper limit of R < 0.41. The column densities, cover- 
ing fractions, and power-law photon index are all consistent at 
the la level with the values in Model IV. The best-fitting pa- 
rameters are given in Table [3] Since the strength of reflection 
correlates with the e-folding energy, we repeated our fits with 
different values of i? cu t, ranging from 100 to 1000 keV. We 
show the upper limits to R as a function of e-folding energy 
in Figure [6] 

We also attempted to fit the spectrum with an ionized 
Compt on reflection model, using the the REFLIONX table 
model (iRoss& Fabianl 120051) (Model VIII). This model re- 
sulted in a similar fit to the spectrum (x 2 =174 for 148 dof), 
and the ionization parameter pegged at its hard limit of £=10 4 
ergs cm s _1 . 

5. FLUORESCENT LINE EMISSION 
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FIG . 6 . — 90% confidence upper limit to the neflection fraction as a function 
of e-folding energy. 
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FIG. 7.— Data/model residuals for Suzaku XIS FI (black) and XIS BI (red) 
spectra in the energy range 4-8 keV, after fitting Model IV to the FI, BI, 
and PIN spectra to the broadband spectrum, but excluding the energy range 
6.2-6.5 keV (rest frame). Prominent residuals at ^6.4 keV are observed, 
indicating the need for a narrow Fe Ko line. No red wing to the Fe Ka 
emission is observed. 
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FIG. 8. — Energy and line width confidence contours for the Fe Ko line 
measured from the Suzaku XIS FI and BI spectra. Shown are the 68%, 90%, 
and 99% contours. 

We next measured the properties of the fluorescent Fe Ka 
emission. To give a qualitative sense of the line profile, we fit- 
ted the continuum-only components of Model IV to the XIS 



FI, XIS BI, and HXD PIN spectra in the energy range 2-70 
keV. We excluded the rest frame Fe Ka bandpass of 6.2-6.5 
keV. The continuum parameters of this fit are consistent at the 
ler level to those found in Model IV. We then extended the 
energy range to include the Fe Ka bandpass. The data/model 
residuals between 4 and 8 keV are shown in Figure [7] There 
are clear residuals at a rest-frame energy of 6.4 keV, corre- 
sponding to neutral Fe Ka emission. 

Ou r best-fitting model (Model IV), which we discussed in 
§ I3. 3 1 includes Fe Ka emission with the following properties: 
its centroid is 6.3851q'q23 keV anc ^ me u PP er limit to the line 
width is <65 eV (both values are quoted at 90% confidence 
for two interesting parameters). This is less than the ^100 eV 
FWHM resolution of the XIS at ~6 keV. Confidence contours 
for the energy and line width are shown in Figure |8] The 
equivalent width of the Fe Ka line is 129 ± 43 eV. The pa- 
rameters of the Fe Ka line are summarized in Table [2] and are 
consistent with those found wi t h Chandra and XM M-Newton 
observations dKraft et alj|2007t lEvans et al.l l2006). The addi- 
tion of a Fe K/3 line is statistically insignificant. 

There is no indication of a velocity broadened component 
of the Fe Ka emission, either in the form of a broad underly- 
ing Gaussian or a relativistically blurred 'diskline'. However, 
given that some radio gal axies may have mod estly broadened, 
ionized Fe Ka lines (e.g., Kataoka et al. 2007), we calculated 
the upper limit to the equivalent width of a Gaussian of frozen 
line width <r=0.5 keV and energy 6.8 keV. We found the 90% 
confidence upper limit to be <65 eV. 

6. SEARCHING FOR VARIABILITY: JOINT FITTING TO SUZAKU, 
SWIFT, CHANDRA, AND XMM-NEWTON SPECTRA 

We searched for inter-epoch spectral variability of 3C 33 by 
simultaneously fitting the Suzaku XIS FI, PI, and PIN, Swift 
BAT, Chandra ACIS (two observations), and XMM-Newton 
EPIC/pn spectra. We adopted our baseline fit (Model IV) for 
our spectral fits. We used this model instead Model HI be- 
cause the normalizations of the soft power laws are likely 
to differ considerably. This is due to the different extrac- 
tion regions used for our Suzaku XIS (2.5'), Chandra (5"), 
and XMM-Newton (35") analyses, which sample dissimilar 
amounts of extended, soft X-ray emission. We restricted the 
energy range of the fits to 0.5-10 keV (Suzaku XIS FI and BI), 
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FIG. 9. — Counts spectra, best fitting model, and residuals for all data. 
Shown are the Suzaku XIS FI (black), XIS BI (red), HXD PIN (green), Swift 
BAT (blue), Chandra OBSID 6910 (light blue), Chandra OBSID 7200 (ma- 
genta), and XMM-Newton pn (orange). The model fit shown is the sum of a 
partially covered power law, a soft, unabsorbed power law, and a Gaussian 
Fe Ka line (Model IV). 
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FIG. 10. — Combined Chandra, XMM-Newton, and Suzaku CCD spectrum of 3C 33 in the energy range 0.5-1 .5 keV. Also plotted are the spectral transitions 
identified in the photoionized plasma of the canonical Seyfert 2 galaxy, NGC 1068. 

15-30 keV (Suzaku PIN), 14-100 keV (Swift BAT), 0.5-7 keV 
(Chandra), and 0.5-8 keV (XMM-Newton). 

For our spectral fitting, we linked the column densities, cov- 
ering fraction and photon index associated with the primary 
power law across all five data sets, and tied the photon index 
of the unabsorbed power law to that of the absorbed com- 
ponent. We allowed the normalizations of the Swift, Chan- 
dra, XMM-Newton, and (combined) Suzaku power laws to 
vary independently, and applied the Suzaku PIN/XIS cross- 
normalization factor of 1.09. Finally, we added an unresolved 
Gaussian Fe Ka line to our model. 

The model provided a good description of all the spectra 
(\ 2 = 292 for 261 dof). All the parameters are shown in Ta- 
ble!?] m summary, the best-fitting parameters of the partially 
covered power law are /Vh=(4.9*q'|) x 10 23 cm~ 2 with a cov- 
ering fraction of 92t 3 %, and photon index T=1.67 ± 0.10. 
The unabsorbed 2-100 keV (extrapolated) luminosities of the 
primary power law are consistent with no inter-epoch vari- 
ability. In Figure [9] we show the counts spectra and model 
residuals for the seven datasets. 

7. SOFT-BAND X-RAY EMISSION 

In order to search for emission lines that might indi- 
cate emissio n from a thermal plasma, as first reported by 
iTorresi et all ([2009), we combined the Chandra, XMM- 
Newton, and Suzaku spectra together using the ISIS spectral 
fitting packagfl We rebinned each data set to a common 
grid, namely the HWHM of the Suzaku XIS CCDs. In Fig- 
ure[l0] we show the fluxed spectrum between 0.5 and 1.5 keV, 
and plot the key photoionization transitions identified in the 
Chandra HET GS spectrum of th e canonical Seyf ert 2 galax- 
ies N GC 1068 (iQgle et al.ll2003h and NGC 2110 (lEvans et alJ 
2006). We interpret these results in § 10 



origins for the soft X-ray emission ( £18.31 ). and finally interpret 
the accretion, reflection and possible jet-production properties 
of 3C 3 3 in the context of other radio-loud and radio-quiet 
AGN ( ^4b . 

8.1. Fe Ka Diagnostics 

The strong Fe Ka line detected with Suzaku allows us 
to place constraints on the location and physical state of 
the fluorescing material. The energy of the line core, 
6.3851 

o'o23 keV, is consistent with Fe fluorescence from neu- 
tral or near-neutral (up to ^Fe XVIII) species. The unre- 
solved Fe Ka fine width and the estimated black hole mass of 
-4 x 10 s M Q (ISmith et al.lll990l : lBettoni et al.ll2003h provide 
a lower limit to the emission radius of 0.1 pc (^2000 R s ), 
using Keplerian arguments. Finally, the equivalent width of 
the Fe Ka line is consistent with transmission through an ab- 
sorbing column of relative iron abundance ^4fo=1-0 and col- 
umn density similar to that measured from o ur spectral fitting 
dGhisellini et al.lll994l : lMiyazaki et al.lll996h . 

There is no evidence for an additional velocity broadened 
component of the Fe Ka line, which rules out the presence 
of relativistically blurred fluorescent emission from the inner- 
most portions of the accretion flow. This may imply that a 
standard, cold Shakura-Sunyaev accretion disk is truncated 
at a large distance from the central supermassive black hole. 
We cannot rule out the presence of a modestly broadened 
(R ~ 100i? s ), ionized Fe Ka line in 3C 33 but nonetheless 
it seems clear that fluorescent emission from the inner disk is 
absent. We will return to this point in ; 



8. INTERPRETATION 

We have shown that the nuclear X-ray spectrum of 3C 33 
is well described by a power law that is transmitted either 
through one or more layers of neutral material, or through a 
partially ionized obscurer. There is no evidence for a Comp- 
tonized reflection continuum above ~10 keV, which means 
we have now resolved the que stion of reflection raised by pre- 
vious Chandra observations (Kra ft et all 120071) . In this sec- 
tion, we use the Fe Ka line to diagno se th e physical location 
and state of the fluorescent emission ( ^8. IK discuss the origin 
of the intrinsic absorption in 3C33 (i j8.2l ). consider various 

9 http://space.mit.edu/cxc/isis/ 



8.2. Origin of Absorption 

We considered in detail several absorption models that 
can account for the nuclear spectrum. In one (Model III), 
two zones of neutral material (a partial coverer of column 
ATh,i~6 x 10 23 cm -2 with a covering fraction of ^80%, plus 
a second column of /Vh.2^2 x 10 23 cm~ 2 ) obscure the pri- 
mary emission. The second column takes the form of a patchy 
absorber if we explicitly choose to associate the soft, unab- 
sorbed, X-ray emission with the primary power law (Model 
IV). Alternatively, we can model the absorber as a single 
zone of modestly ionized material (Model V), which transmits 
more of the low-energy continuum emission than a neutral ob- 
scurer and thus reproduces the spectral complexity in the 2-10 
keV spectrum. The relatively mild ionization of the warm ab- 
sorber is also consistent with the neutral or near-neutral Fe Ka 
emission observed. Again, the warm absorber is patchy if we 
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choose to associate the soft X-ray continuum with the same 
physical process as the primary power law. 

Both models have distinct advantages over simple param- 
eterizations of the nuclear continuum: they adequately ex- 
plain the low-energy spectrum and reconcile the photon index 
of the primary power law with canonical values found in ra- 
dio galaxies. However, we are unable to distinguish between 
these models due to the lack of photon statistics at lower en- 
ergies, where discrete ionized absorption features would be 
observed in the case of the warm absorber. 

If we choose to associate the unresolved fluorescent Fe Ka 
emission with the circumnuclear absorber then, as calculated 
above, the lower limit to the radius of the obscurer is ~0.1 pc. 
We can also calculate the upper limit to the radius, r, in the 
case of the warm absorber model, since the ionizing lumi- 
nosity, L = (nr 2 , where £ is the ionization parameter and 
n is the electron number density. Assuming that the thick- 
ness of the absorber is less than the distance to the absorber 
from the source of ionization, i.e., Ar/r < 1, we can estab- 
lish the inequality r < L/Nn£. We find R 55 20 pc, which 
is consistent with the parsec-scale absorber invoked in AGN 
unification models. 

8.3. Origin of Soft X-ray Emission 

There is clear evidence of a 'soft excess' in the X-ray spec- 
trum of 3C33. We consider several possibilities for its ori- 
gin here. First, the emission could be the result of patchy 
circumnuclear absorption (either neutral or partially ionized), 
which allows a small fraction of continuum emission to be 
transmitted free from obscuration. Such a model is consistent 
with observations of numerous AGN, including NGC2110 
(Evans et al. 2007). The lack of variability of the source on 
timescales of years (SQ suggests that such a patchy absorber 
would subtend a large solid angle at the source of the X-ray 
emission. 

Second, the soft emission could represent the scattering of 
the primary (heavily obscured) power-law emission in to the 
line of sight, as is pred icted by unified models of AGN (e.g., 
IUrrv & PadovanIl [l995h . Here, the emission is scattered by 
free electrons in the vicinity of the broad-line region and may 
be associated with the (possibly photoionized) circumnuclear 
regions of the AGN. The 'efficiency' of scattering is simply 
given by the ratio of the normalizations of the primary and 
soft power laws — for 3C33, we find the efficiency to be 
~2%. This number is consisten t with the typical r ange found 
in Seyfert galaxies of 1-10% (iTurner et al.l Il997h . Further 
support for this model comes from the spatially extended soft 
X-ray and [O III] on scales of ^5 kpc emis sion seen in Chan - 
dra and HST observations of the source (lKraftetal.ll2007h . 
The emission is elongated along the axis of the large-scale 
radio emission in the source, and therefore likely along the 
axis of the circumnuclear obscurer described by unification 
models. In this case, we might expect a contribution from 
photoionized emission lines in the soft X- r ay spe ctrum in the 
case of 3C 33, as argued by Torr esi et all (120091) . Our com- 
bined soft X-ray spectrum (Fig. ITQb shows possible hints of 
He-like species of O VII and Mg XI, but there are insufficient 
photon statistics to determine if we are witnessing photoion- 
ized gas. 

A third and final possibility is that the soft X-ray emis- 
sion is not associated with the primary X-ray emission re- 
lated to the accretion flow, and instead originates in an unre- 
solved parsec-scale jet. Given that 3C 33 is a strongly radio- 
loud AGN, a contribution to the X-ray emission from the 
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FIG. 1 1 . — 1-keV X-ray luminosity of the soft X-ray emission component 
as a function of 5-GHz radio core luminosity in all observed 3CRR sources at 
z < 0.5 lEvans et al. 2006; Hardcastle et al. 2006). Open circles are LERG, 
filled circles NLRG, open stars BLRG, and filled stars quasars. Large sur- 
rounding circles indicate that a source is an FRI. 3C 33 is marked by a box. 
Dotted lines show the regression line to all data and its 1<t confidence range. 
The similarity of 3C 33 to the other sources suggests that its soft X-ray com- 
ponent originates in an unresolved jet, although we cannot rule out other 
interpretations. 

jet is not unexpected, and has been argued for extensively 
in o ther radio galaxies by, e. g., [Hardcastle & WorralU d 19991) 
and iBalmaverd eet al.l (120061) . The key to these authors' in- 
terpretations is the strong observed correlation between the 
X-ray and radio core luminosities in large samples of radio- 
loud AGN, which they argue supports a common origin of the 
two components at the base of a jet. In FigureQTl we show the 
1-keV and 5-G Hz luminosity densities of all observed z<0.5 
radio galaxies (Hardcastl e et al J 120061) . together with 3C33. 
This figure shows that 3C 33 has a similar measured X-ray 
luminosity to that predicted from the lHardcastle et aTl (120061) 
relation. 

In summary, we cannot distinguish between these three 
physical scenarios. Future observations with the gratings in- 
struments on board Chandra and XMM-Newton would allow 
us to determine which is correct. 

8.4. Accretion and Reflection in 3C33 

Our results demonstrate that 3C 33 shows no signs of 
Compton reflection from neutral material in the inner re- 
gions of an accretion disk: there is no reflection hump at en- 
ergies >10 keV, and no evidence for relativistically broad- 
ened Fe Ka emission (Model VII). Indeed, the narrow Fe 
Ka line is likely to originate in a much farther region, at 
least 2,000 R s from the black hole. Further, as 3C 33 is a 
narrow-line radio galaxy, meaning we can spectrally sepa- 
rate the X-ray emission associated with the accretion flow 
and jet, we rule out the idea that X-ray emission from the 
jet dilutes the reflection spectrum of the source at high ener- 
gies. The lack of a Compton hump in 3C 33 is consistent with 
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growing numbers of radio-loud AGN observed with Suzaku to 
show systematical ly weaker reflection than their radio-quiet 
counterparts (e .g., Larsson et al. 2008; Samb runa et alJl2009l : 
iKataoka eT al. 2007). There are two principal models that can 
explain this dichotomy, which basically stem from differences 
in either the accretion flow geometry, or the black-hole spin 
configuration. We examine both in the context of 3C 33 be- 
low. 

8.4.1. Ionized Accretion Models 

iBallantyne et al.l (l2002h interpret the weak reflection fea- 
tures observed in radio-loud AGN such as 3C 33 to be the 
result of an ionized inner accretion flow. Since the reflected 
spectrum of an incident power law on to highly ionized mate- 
rial is itself a power law, they demonstrate that the reflection 
fraction, R, can remain ar bitrarily high. Our attempt to fit an 
ionized reflection model (Ross & Fabian] 120051) to 3C33 led 
to the ionization parameter pegging at its hard limit of £=10 4 
ergs cm s _1 (Model VIII). For 3C33, therefore, we cannot 
distinguish between a very highly ionized reflector, or simply 
an absent inner disk. We further note that the lack of ionized 
Fe Ka emission is consistent with this hypothesis. 

8.4.2. Retrograde Black-Hole Spin 

The X-ray spectrum of 3C 33 is also consis tent with the pic- 
ture of Garofalo, Ev ans. & Sambrurial ((2009). They argue the 
jet power of high-excitation radio galaxies such as 3C 33 is 
the result of retrograde black-hole spin with respect to the 
accreting material. This configuration results in a larger in- 
nermost circular stable orbit than for a prograde black hole. 
The weak or absent signatures of reflection in 3C 33, then, 
are simply a consequence of the large inner dis k radius. Fur- 
thermore, Garofalo, Evans, & Sambruna d2009l) show that the 
retrograde spin configuration produ ces both strong Blandford- 
Pay ne and Blandford-Znajek j ets dBlandford & Paynd[l982l 
and Blandford & Znajek 1977, respectively). They point out 
that this highly relativistic motion away from the inner ac- 
cretion flow can easily wash out reflection fea tures, as first 
demon strated by Reynolds & Fabian] (1 1997b and IB eloborodovl 
(1999). Our X -ray observations of 3C33 support this sce- 
nario. 

In summary, we have presented two, equally plausible, 
models that can explain the paucity of a Compton reflection 
hump in 3C 33. Future, high-sensitivity observations of 3C 33 
with IXO will be needed to distinguish between them. 

9. CONCLUSIONS 

We have presented results from a new 100-ks Suzaku ob- 
servation of the radio galaxy 3C33. The Suzaku data have 
allowed us to break degeneracies in the reflection and ab- 
sorption properties of 3C 33 that were present in our previous 
Chandra and XMM-Newton observations. We have shown the 
following: 

1. The nuclear X-ray spectrum of 3C 33 is well described 
by a power law of photon index ~1.7 that is absorbed 
either through one or more layers of pc-scale neutral 
material, or through a partially ionized pc-scale ob- 
scurer. Both models have distinct advantages over sim- 
ple parameterizations of the nuclear continuum: they 
adequately explain the low-energy spectrum and recon- 
cile the photon index of the primary power law with 
canonical values found in radio galaxies. 



2. There is evidence for a 'soft excess' in the X-ray spec- 
trum of 3C 33. Its origin is consistent with several pos- 
sibilities, including patchy circumnuclear absorption, 
photoionized gas, and an unresolved parsec-scale jet. 

3. 3C 33 shows no signs of Compton reflection from neu- 
tral material in the inner regions of an accretion disk: 
the upper limit to the reflection fraction for neutral ma- 
terial is R < 0.41 for an e-folding energy of 1 GeV 

4. The observed Fe Ka line originates in neutral or near- 
neutral material at least 2,000 R s from the black hole. 
We cannot rule out the presence of a modestly broad- 
ened (R ~ 100i? s ), ionized Fe Ka line in 3C 33 but 
nonetheless it seems clear that fluorescent emission 
from the inner disk is absent. 

5. The absence of a Compton reflection hump in 3C33 
is consistent with eith er an ionized accretion flow 
(Ball antvne et al.ll2002l). or with retrograd e black-hole 
spin (Garofalo, Evans, & Sambruna 2009). Future ob- 
servations with IXO will allow us to distinguish be- 
tween these two models. 
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TABLE 1 
Observation Log 



Observatory 


Date of observation 


Obs. ID 


Instrument 


Screened exposure time 


Net nuclear count rate (s ) 


Suzaku 


2007 Dec 26 


702059010 


XIS0 


127 ks 


(5.47 ±0.14) x 10~ 2 




2007 Dec 26 




XIS1 


127 ks 


(5.07 ±0.14) x 10~ 2 




2007 Dec 26 




XIS3 


127 ks 


(5.17 ±0.08) x 10~ 2 




2007 Dec 26 




HXD/PIN 


99 ks 


(5.22 ±0.26) x 10" 2 


Swift 


2004 Dec 15 - 2006 Nov 01 


SWIFTJ0109.0+1320 


BAT 


2.8 Ms (effective) 


(6.31 ±0.89) x 10~ 4 


Chandra 


2005 Nov 08 


6910 


ACIS-S 


20 ks 


(6.38 ±0.16) x 10~ 2 


Chandra 


2005 Nov 12 


7200 


ACIS-S 


20 ks 


(6.55 ±0.18) x 10" 2 


XMM-Newton 


2004 Jan 21 


0203280301 


EPIC/pn 


6.3 ks 


(1.37 ±0.05) x 10" 1 



TABLE 2 
Fe Kq properties from Model IV 



Parameter 


Value 




Rest-frame energy 


6.385±g;gl| keV 




Line width 


<65 eV [XIS FWHM 


= lOOeV] 


Intensity 


(2.40 ±0.04) x 10~ 5 


photons cm~ 2 s _1 


Equivalent width 


129±43 eV 





TABLE 3 

Continuum Spectral Parameters for Joint Model Fits to Suzaku and Swift Spectrum. Luminosities quoted are unabsorbed and extrapolated to 2-100 keV. 











Absorber 


Covering 






Reflector 


-^(2-lOOkeV) 
(Power Law) 




Overview 


Model 


Description of spectrum 


N n (cm" 2 ) 


log(0 


fraction 


r 


H 


logCO 


(ergs s _1 ) 


X 2 /dof 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Single 


I 


Abs(PL( 1 )+Gauss)+PL(2) 


(3.7 ±0.20) x 10 23 




100% 


ri=1.39±0.05 






(3.9+ lJ ;e) x 10 44 


242/152 


absorption 












r 2 =ri 






(1.7 ±0.1) x 10 43 




Multiple 


II 


PCAbs(PL+Gauss) 


(3.7 ±0.20) x 10 23 




96±1% 


1.39 ±0.05 






(4.1+^) x 10 44 


242/1 5-^ 


absorption 


III 


PC Abs( 1 ) X PCAbs(2) X (PL+Gauss) 


AT H ,l=(6.2t? ; 8) x 10 23 
AT H 2=(l-6 ± 0.4) x 10 23 




84+? % 
91+| % 


1.71 ± 0.09 






(4.5+il) x 1044 


176/15^. 
p 




IV 


PCAbs( 1 ) X Abs(2) X (PL( 1 )+Gauss)+PL(2) 


JVh!i=(6.01;;») x 10 23 




85+4 % 


ri=1.71 ± 0.09 






(4A+H) x 10 44 


176/15^ 








Af H ,2=(l-7±0.4) x 10 23 




100% 










1 

(a 














r 2 =ri 






(6.5 ±0.4) x 10 42 


Warm 


V 


PCWarmAbs(PL( 1 )+Gauss) 


(5.2^; 3 ) x 10 23 


1 ' ol -0.20 


98±1% 


ri=1.67±0.07 






(4.2+^;'() x 10 44 


184/15^- 


absorption 






(5.2±0; 3 ) x 10 23 


1 40+017 










(4.1±g;|) x 10 44 
(6.9 ±0.4) x 10 42 


S' 




VI 


WarmAbs(PL( 1 )+Gauss)+PL(2) 


100% 


ri=1.69±0.07 

r 2 =ri 






184/15? 
O 




Neutral 


VII 


PC Abs( 1 ) X Abs(2) X (RefPL( 1 )+Gauss)+PL(2) 


N H i=(5.7± 1.0) x 10 23 




86±6% 


Ti=1.75 ± 0.08 


<0.41 




(4.0+ 1 / 4 ) x 10 44 
(5.8 ± 0.4) x 10 42 


175/142) 


reflection 






Nu,2=(^° i) x io 23 




100% 


r 2 =ri 








Ionized 


VIII 


PCAbs(l) X Abs(2) X (PL(l)+Reflion(2)+Gauss)+PL(3) 


JVh,i=(6.6 ± 0.2) x 10 23 




83±1% 


1 1-1-O»_o.02 






(1.9 ±0.2) x 10 44 


174/148 


reflection 






JV H ,a=(2.0±0.1) x 10 23 




100% 


r 2 =r! 
r 3 =ri 




4 (pegged) 


(2.5 ±0.1) x 10 44 
(9.6 ± 0.5) x 10 42 





Col. (1): Model number. Col. (2): Description of spectrum (Abs=Neutral absorption, PCAbs=Partially covering neutral absorption, WarmAbs=Warm Absorption (£ is the ionization parameter and is quoted in units 
of ergs cm s ), PCWarmAbs=Partially covering warm absorption, PL=Power Law, RefPL=Sum of direct and reflected power-law continuum from neutral material, Reflion= Ross & Fabian 1 2005) ionized reflection 
continuum, Gauss=Redshifted Gaussian line. Col. (3): Intrinsic neutral hydrogen column density. Gal actic absorption has also been applied. Col. (4): Ionization parameter of absorber. Col. (5): Covering fraction. 
Col. (6): Power-law photon index. Col. (7): Reflection fraction. Col. (8): Ionization parameter for Ross & Fabian (2005) ionized reflection model. Col. (9): 2-100 keV (extrapolated) unabsorbed luminosity of 
primary power law. Col. (10): Value of \ 2 and degrees of freedom. 
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TABLE 4 

Best-fitting parameters for joint fit to Suzaku, Swift, Chandra, and XMM-Newton Spectra. Luminosities quoted are 

unabsorbed and extrapolated to 2-100 kev. 



Component 



N H (cm- 2 ) 



Covering 
fraction F 



(9.8t*; 4 ) x 10 22 



-^(2-lOOkeV) 

(Power Law) 
(ergs s -1 ) 



ichandra,6910 = (■i^-li's) X 1Q 44 



Partially covering 
primary power law 



(4.9+|l) x 10 2 



100% 



1.67 ±0.10 



ichandra,7200 = (3.5^q;°) X 10 4 

ixMM = (3.3^ 'g) x 10 44 



MM = (3.3 j 

-^Suzaku 

Lswift = (4.1 



(3.8tJ; 2 

-1.9 \ 
-1.2) 



x 10* 
x 10 44 



Soft power law 



100% =r x L Chandrai6910 = (2.1 ± 0.3) x 10 42 

Lchandra,7200 = (2.5 ± 0.4) X 10 44 

ixMM = (3.5 ± 0.5) x 10 44 
Lsuzaku = (7.0 ± 0.5) x 10 44 



